function u = Lax_Friedrichs(u,N,mu,num)
for i = 1:num
    u = Lax_Friedrichs_onestep(u,N,mu);
end